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A Mathematical Model for the Behavior of 
Individuals in a Social Field 



Abstract 

Related to an idea of Lewin, a mathematical model for behavioral changes under the 
influence of a social field is developed. The social field reflects public opinion, social 
norms and trends. It is not only given by external factors (the environment) but also 
by the interactions of individuals. Two important kinds of interaction processes are 
distinguished: Imitative and avoidance processes. Variations of individual behavior 
are taken into account by "diffusion coefficients" . 

Key words: Behavioral model, field theory, diffusion model, decision theory, pair 
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A Mathematical Model for the Behavior of 
Individuals in a Social Field 

1 Introduction 

Many models have been developed for behavioral changes, but only a few are formulated 
in terms of mathematical relations. For example, 

• game theory (von Neumann and Morgenstern, 1944), based on the concept of 
success of meeting strategies, is used for the description of cooperation and compe- 
tition processes among individuals, 

• decision theories (Domencich and McFadden, 1975, Ortuzar and Willumsen, 
1990), assuming the maximization of utility, successfully model the choice behavior 
among several alternatives, 

• diffusion models (Coleman, 1964, Bartholomew, 1967, Granovetter, 1983, 
Kennedy, 1983) mathematically describe the spread of behaviors or opinions, ru- 
mors, innovations, etc. 

All these models are related to a more general behavioral model discussed in the fol- 
lowing. This model is based on BoLTZMANN-/zA;e equations and includes spontaneous (or 
externally induced) behavioral changes and behavioral changes by pair interactions of 
individuals (sect. |2|). These changes are described by transition rates. They reflect the 
results of mental and psychical processes, which could be simulated with the help of Os- 
good and Tannenbaum's (1955) congruity principle, Heider's (1946) balance theory or 
Festinger's (1957) dissonance theory. However, it is sufficient for our model to deter- 
mine the transition rates empirically (sect. |5|). The ansatz used for the transition rates 
distinguishes imitative and avoidance processes, and assumes utility maximization by a 
variant of the multinomial logit model (Domencich and McFadden, 1975, Ortuzar 
and Willumsen, 1990) (sect. gjg ). 

In section || a consequent mathematical formulation related to an idea of Lewin (1951) is 
developed, according to which the behavior of individuals is guided by a social field. This 
formulation is achieved by a second order Taylor approximation of the Boltzmann- 
like equations leading to diffusion equations. Because of their relation with the Boltz- 
mann equation (Boltzmann, 1964) and the Fokker-Planck equation (Fokker, 1914, 
Planck, 1917) they will be called the Boltzmann-Fokker-Planck equations. Accord- 
ing to these new equations the most probable behavioral change is given by a vectorial 
quantity that can be interpreted as social force (sect. |3.1| ). The social force results from ex- 
ternal influences (the environment) as well as from individual interactions. In special cases 
the social force is the derivative (gradient) of a potential. This potential reflects public 
opinion, social norms and trends, and will be called the social field. By diffusion coeffi- 
cients individual variation of the behavior (the "freedom of will") is taken into account. 
In section |] representative cases are illustrated by computer simulations. 
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The Boltzmann-Fokker-Planck model for the behavior of individuals under the in- 
fluence of a social field shows some analogies with the physical model for the behavior of 
electrons in an electric field (e.g. of an atomic nucleus) (Helbing, 1992a,c). In particular, 
individuals and electrons influence the concrete form of the effective social, respectively, 
electric field. However, the behavior of electrons is governed by a different equation: the 
Schrodinger equation (Schrodinger, 1926, Davydov, 1976). 



2 The BOLTZMANN-like behavioral model 

Let us consider a population consisting of a great number N ^> 1 of individuals. Concern- 
ing a special topic of interest, these individuals show a behavior x out of several possible 
behaviors in the set Q. 

Due to "freedom of the will" one cannot expect a deterministic theory for the temporal 
change dx/dt of the individual behavior x(t) to be realistic. However, one can construct 
a model for the change of the probability distribution P(x,t) of behaviors x(t) within 
the given population (P(x,t) > 0, P(x,t) = 1). A theory of this kind is, of course, 

stochastic. In order to take into account several types a of behavior, we may distinguish 

A 

A subpopulations a consisting of N a ^> 1 individuals (Yj N a = N). Then, the following 

a=l 

relation holds: 

P{x,t)=i:JjfPa{x,t). (1) 

Our goal is now to find a suitable equation for the probability distribution P a (x,t) of 
behaviors within subpopulation a (P a (x,i) > 0, P a (x,t) = 1). If we neglect memory 

effects (cf. sect, p.lp , the desired equation is of the form 

—P a (x, t) = inflow into x — outflow from x . (2) 

Whereas the inflow into x is given as the sum over all absolute transition rates describing 
changes from an arbitrary behavior x' to x, the outflow from x is given as the sum over 
all absolute transition rates describing changes from x to another behavior x'. Since the 
absolute transition rate of changes from x to x' is the product w a (x'\x; t)P a (x, t) of the 
relative transition rate w a (x'\x; t) for a change to behavior x' given x, and the probability 
P a (x,t) of behavior x, we arrive at the explicit equation 

j t P a (x,t)= £ [w\x\x'-t)P a {x\t)-w\x'\x-t)P a {x,t)\. (3) 

x 1 en 

w a (x'\x; t) has the meaning of a transition probablility from x to x' per unit time and 
takes into account the behavioral variations between the individuals (occurring even 
within the same type a of behavior!). 

In the following we have to specify the relative transition rates w a (x'\x; t), which will turn 
out to be effective transition rates. If we restrict the model to spontaneous (or externally 
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induced) behavioral changes and behavioral changes due to pair interactions, we have 
(Helbing, 1992a,c): 



u>°(se'|a:;t) := w a (x'\x;t) E N b w ab (x',y'\x,y;t)P b (y,t) 

6=1 yen y' en 



(4) 



w a (x'\x;t) describes the rate of spontaneous (resp. externally induced) transitions from 
x to x 1 for individuals of subpopulation a. w ab (x', y'\x, y\ t) is the transition rate for two 
individuals of types a and b to change their behaviors from x and y to x' and y' , respec- 
tively, due to pair interactions. The total frequency of these interactions is proportional 
to the probability P b (y,t) of behavior y within subpopulation b and the number N b of 
individuals of type b. We have to sum up over b, y, and y' since all specifications of 
these variables are effectively connected with transitions from x to x' of individuals of 
subpopulation a. 



Inserting 
1992a,c) 



into (|3p, we now obtain the socalled Boltzmann- like equations (Helbing, 
d 



dt 



(5a) 



with 



P a (x,t) = [w a {x\x';t)P a (x',t) - w a {x'\x;t)P a (x,t) 

x'en 

A 

+ w ab (x,y'\x',y;t)P b (y,t)P a (x>,t) 

6=i x'en yen y'en 

A 

~ E E E E w ab(x', y'\x, y; t)P b (y, t)P a {x, t) (5b) 
6=i x'en yen y'en 

w ab (x',y'\x,y;t) := N b w ab (x', y'\x, y; t) . (6) 

Obviously, (§d) depends nonlinearly (quadratically) on the probability distributions 
P a (x,t) (resp. P b (y,t)) which is due to the pair interactions. 

The BoLTZMANN-like equations originally had been developed for the description of the 
kinetics of gases (Boltzmann, 1964). However, they have also been applied to atti- 
tude formation (Helbing, 1992b, c) and the avoidance behavior of pedestrians (Helbing, 
1992c,d). 

It is possible to generalize the model to simultaneous interactions of an arbitrary number 
of individuals (i.e., higher order interactions) (Helbing, 1992a,c). However, in most cases 
behavioral changes are dominated by pair interactions (dyadic interactions). Many of the 
phenomena occurring in social interaction processes can already be understood in terms 
of pair interactions. 



2.1 The form of the transition rates 



For models of behavioral changes the following special form of the effective transition 
rates (f|) has been found to be suitable (Helbing, 1992b, c,e): 

A 

w a (x'\x; t) := u a {t)R a (x'\x; t) + £ u ab (t) \fl b (t)P b (x\ t) + f 2 ab (t)P b (x, tj\ R a (x'\x; t) . (7) 



6=1 



3 



Here, 



• u a (t) is a measure of the rate of spontaneous (or externally induced) behavioral 
changes within subpopulation a. 

• R a (x'\x; t) [resp. R a (x'\x; t)) is the readiness of an individual of subpopulation a to 
change behavior from x to x' spontaneously [resp. in pair interactions]. 

• v a b(t) = NbVabit) is the interaction rate of an individual of subpopulation a with 
individuals of subpopulation b. 

• fa b (t) is a measure for the frequency of imitative processes 

x', x' < — x, x [x 7^ x') , (8) 

where an individual of subpopulation a takes over the behavior x' of an individual 
of subpopulation b. The total frequency of imitative processes is proportional to the 
probablility Pb(x' ,t) of behavior x' within subpopulation b. 

• fab(t) is a measure for the frequency of avoidance processes 

x', x < — x,x (x 7^ x 1 ) , (9) 

where an individual of subpopulation a changes the behavior x to another behavior 
x' if meeting an individual of subpopulation b with the same behavior (defiant 
behavior, snob effect). The total frequency of avoidance processes is proportional to 
the probablility Pb(x,t) of behavior x within subpopulation b. 

A more detailled discussion of the different kinds of interaction processes and of ansatz 
(0) is given in publications of Helbing ( 1992b, c,e). 

For R a (x'\x; t) we take the quite general form 

JJ a {x' ,t)-U a (x,t) 

Da{x , x . t) (10a) 

with 

D a (x',x;t) = D a (x,x';t) > 

(cf. Weidlich and Haag, 1988, Helbing, 1992c). Then, the readiness R a (x'\x; t) for an 
individual of subpopulation a to change behavior from x to x' will be greater, 

• the greater the difference in the utilities U a (.,t) of behaviors x' and x, 

• the smaller the incompatibility ("distance") D a (x',x;t) between the behaviors x 
and x'. 



Similar to ( |10a| ) we use 

e u a (x' ,t)-U a {x,t) 

= < 10b > 

4 



and, therefore, allow the utility function U a (x,t) for spontaneous (or externally induced) 
behavioral changes to differ from the utility function U a (x,t) for behavioral changes in 
pair interactions. Ansatz ( |T0D is related to the multinomial logit model (Domencich and 
McFadden, 1975, Ortuzar and Willumsen, 1990). It assumes utility maximization 
with incomplete information about the exact utility of a behavioral change from x to x', 
which is, therefore, estimated and stochastically varying (cf. Helbing, 1992c). 

Computer simulations of the BoLTZMANN-like equations (|3|), (0), ( |i"0|) are discussed and 
illustrated in Helbing (1992b,c,e) (cf. also sect. |). 



2.2 Special fields of application in the social sciences 

The BoLTZMANN-like equations (|3|), (|7|) include a variety of special cases, which have 
become very important in the social sciences: 

• The logistic equation (Pearl, 1924, Verhulst, 1845) describes limited growth 
processes. Let us consider the situation of two behaviors x = x G {1,2} (i.e., 
P a (l,t) = 1 — P a {2,t)) and one subpopulation (A = 1). x = 2 may, for example, 
have the meaning to apply a certain strategy, and x = 1 not to do so. If only 
imitative processes 

2,2 <— 1,2 (11) 
and processes of spontaneous replacement 

1 < — 2 (12) 

are considered, one arrives at the logistic equation 

jPi{%t) = -z/ 1 (t)i? 1 (l|2;t)P 1 (2,t) + z/ 11 (t)/ 1 1 1 (t) J R 1 (2|l;t)(l-Pi(2,t))P 1 (2,t) 
= A(t)P 1 (2,t)(B(t)-P 1 (2,t)). (13) 



The gravity model (Zipf, 1946, Ravenstein, 1876) describes processes of exchange 
between different places x. Its dynamical version results for R a (x'\x] t) = 0, fl b {t) = 
1, f 2 ab {t) = 0, and A — I: 



e U(x,t)-U(x',t) e U(x',t)-U(x,t) 



D(x,x') D(x'x) 



P(x,t)P(x',t). (14) 



Here, we have dropped the index a because of a = 1. P(x,t) is the probability of 
being at place x. The absolute rate of exchange from x to x' is proportional to the 
probabilities P(x, t) and P(x', t) at the places x and x'. D(x, x') is often chosen as a 
function of the metric distance \\x — x'\\ between x and x'\ D(x,x') = D(\\x — x'\\). 

The behavioral model of Weidlich and Haag (1983, 1988, Weidlich, 1991, 1994) 
is based on spontaneous transitions. We obtain this model for /„ 6 (i) = ee f^ b {t) 
and 

A 

U a (x, t) := 6 a (x, t) + Y, K ab P b (x, t) . (15) 
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Because of the dependence of the utilities U a (x,t) on the behavioral distributions 
Pi)(x,t) the model assumes indirect interactions, which are, for example, mediated 
by the newspapers, TV or radio. S a (x,t) is the preference of subpopulation a for 
behavior coupling parameters describing the influence of the behaviorial 

distribution within subpopulation b on the behavior of subpoplation a. For > 0, 
K a b reflects the social pressure of behavioral majorities. 

The game dynamical equations (Hofbauer and Sigmund, 1988, Schuster et. al, 
1981, Helbing, 1992c,e, 1993) result for f^(t) = 5 ab , fl b (t) = 0, and 



R a (x'\x; t) := max \E a (x', t) - E a (x, t), 0J , (16) 

where 

. f 1 if a = b , , , f x if x > y „„ s 

dab := < r. . £ , , and maxlx, y) := < . r (17) 

I if a ^ b v ' y; \ y liy > x. v; 

For a detailled interpretation of these relations see Helbing ( 1992c, e). 
The explicit form of the game dynamical equations is 

^rP a (x,t) = J2[Mx\x';t)P a (x',t)-w a (x'\x;t)P a (x,t)] (18a) 
+ u aa (t)P a (x,t)[E a (x,t) - (E a )] . (18b) 

Whereas flTS|a) again describes spontaneous behavioral changes ( "mutations", inno- 
vations), (]18|b) reflects competition processes leading to a "selection" of behaviors 
with a success E a (x, t) that exceeds the average success 

(E a ) := Y, E a (x',t)P a (x',t). (19) 

x'en 

The success E a (x,t) is connected with the socalled payoff matrices A ab = 
(A ab (x, y)) by 

A 

E a {x, t) := A a (x) + E E A ab{*, y)Pb(y, t) (20) 

6=1 yen 

(Helbing, 1992c,e). A a (x) means the success of behavior x with respect to the 
environment. 

Since the game dynamical equations (0) agree with the selection mutation equa- 
tions (Hofbauer and Sigmund, 1988) they are not only a powerful tool in social 
sciences and economy (Axelrod, 1984, von Neumann and Morgenstern, 1944, 
Schuster et. al, 1981, Helbing, 1992c,e, 1993), but also in evolutionary biology 
(Fisher, 1930, Eigen, 1971, Eigen and Schuster, 1979, Feistel and Ebeling, 
1989). 
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3 The Boltzmann-Fokker-Planck equations 



We shall now assume the set Q of possible behaviors forms a continuous space. The n 
dimensions of this space correspond to different characteristic aspects of the considered 
behaviors. In the continuous formulation, the sums in (^), (|J) have to be replaced by 
integrals: 



^P.(x,t) 



d n x' 



w a (x\x'; t)P a (x', t) - w a (x'\x; t)P a (x, t) 



(21a) 



w a {x'\x;t) :=w a (x'\x;t)+J2 J d n y J d n y' N b w ab (x' ,y'\x,y;t)P b {y,t) . (21b) 

A reformulation of the BoLTZMANN-like equations (|21~1) via a second order Taylor ap- 
proximation (Kramers-Moyal expansion (Kramers, 1940, Moyal, 1949)) leads to 
diffusion equations (Helbing, 1992a, c): 











1 



d d 



^ «) = - E ^ [«■-(*. «)J + 5 E aii dx 
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Q aij (x,t)P a (x,t) (22a) 



with the effective drift coefficients 

K ai (x, t) := J d n x' (x[ - Xi)w a (x'\x; t) 
h 

and the effective diffusion coefficients^ 



(22b) 



(22c) 



Because of their relation with the Boltzmann equation (Boltzmann, 1964) and the 
Fokker-Planck equation (Fokker, 1914, Planck, 1917) equations (fj|) will be called 
the Boltzmann-Fokker-Planck equations in the following (cf. Helbing 1992a,c). In 
the Boltzmann-Fokker-Planck equations the drift coefficients K ai (x,t) govern the 
systematic change ("drift", motion) of the distribution P a (x,t), whereas the diffusion 
coefficients Q a ij(x,t) describe the spread of the distribution P a (x,t) due to fluctuations 
resulting from the individual variation of behavioral changes. 

For ansatz (|7|), the effective drift and diffusion coefficients can be split into contributions 
due to spontaneous (or externally induced) transitions (k = 0), imitative processes (k = 
1), and avoidance processes (k = 2): 



K at (x, t) = £ K k Jx, t) , Q aij (x, t) = £ Q k aij (x, t) 



(23a) 



fc=0 



fc=0 



1 In a paper of Helbing (1992a) the expression for Q a ij{x, t) contains additional terms due to another 
derivation of (p^). However, they make no contributions, since they result in vanishing surface integrals 
(Helbing, 1992c). 
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where 



and 



K 2 ai (x,t) 



v a (t) J d n x' (x[ - Xi)R a (x'\x;t) , 

£ / dV (a/ { - x^R^x'lx; t)P b (x', t) , 

6=1 J 

A 

T.^{t)f ab {t) Jd n x' (xJ-zOi2VI*; WOM) 



<3°j,-(a;, t) := z/„(t) / cfV (x- - x i )(x , j - Xj)R a (x'\x; t) 



(23b) 



A 

Qlijfat) ■= J2^ b (t)f ab (t) <rj(x' i -x i )tf j -x j )ir(x , \x',t)i\(x , ,t), 

6=1 J 



The behavioral changes induced by the environment are included in K^(x,t) and 



3.1 Social force and social field 

The Boltzmann-Fokker-Planck equations ( p2|) are equivalent to the stochastic equa- 
tions (Langevin equations, 1908) 



dx- - 
— - = F ai (x, t) + £ G aij (x, t)^(t) 
at j=1 



with 



and 



1 n 

K ai {x,t) = F ai {x,t) + - £ 

j,fc=l L 



d_ 

dx k 



G ajk (x,t) 



Q aij (x,t) = J2 G aik(x,t)G akj (x,t) 

k=l 



(24a) 



(24b) 



(24c) 



(cf. Stratonovich, 1963, Helbing, 1992c). For an individual of subpopulation a the 
vector £ a (x,t) with the components 



(25) 



describes the contribution to the change of behavior a? that is caused by behavioral fluctu- 
ations £(£) (which are assumed to be delta-correlated and GAUSSian (Helbing, 1992c)). 
Since the diffusion coefficients Q a ij(x,t) and the coefficients G a ij{x,t) are usually small 
quantities, we have F ai (x,t) « K ai (x,t) (cf. ( |24b| ) ) , and ( |24a|) can be put into the form 



— ~ K a (x,t) + fluctuations. 



(26) 
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Whereas the fluctuation term describes individual behavioral variations, the vectorial 
quantity 

/ K al (x,t) \ 

K a (x,t): = : (27) 

\ K an (x,t) J 

drives the systematic change of the behavior x(t) of individuals of subpopulation a. There- 
fore, it is justified to denote K a (x, t) as social force acting on individuals of subpopulation 
a. With that we have attained a very intuitive formulation of social processes, according 
to which behavioral changes are caused by social forces. 

On the one hand, social forces influence the behavior of the individuals, but on the other 
hand, due to interactions, the behavior of the individuals also influences the social forces 
via the behavioral distributions P a (x,t) (cf. (pT| b) , (|22b| )). That means, K a (x,t) is a 
function of the social processes within the given population. 



Under the integrability conditions 

d 



d 



K ai (x, t) = —K aj (x, t) for all i, j 

OXi 



there exists a time- dependent potential 

X n X 

V a (x,t) :=- j dx'-K a (x',t) = -J2J dx'iKaiix 1 ', t) 



i=i 



so that the social force is given by its derivative (by its gradient V): 





K a (x,t) = -VV a (x,t) 



i.e., K ai (x,t) 



V a (x,t) 



(2f 



(29) 



(30) 



The potential V a (x,t) can be understood as social field. It reflects the social influences 
and interactions relevant for behavioral changes: the public opinion, trends, social norms, 
etc. 



3.2 Discussion of the concept of force 

Clearly, the social force is not a force obeying the NEWTONian (1687) laws of mechanics 
(cf. Greenwood, 1988). Instead, the social force K a (x,t) is a vectorial quantity with 
the following properties: 

• K a (x,t) drives the temporal change dx/dt of another vectorial quantity: the be- 
havior x(t) of an individual of subpopulation a. 

• The component 

K ab (x,t) := u ab (t) J d n x' (x' - x)[/i(t)n(aj , ,t) + f 2 ab (t)P b (x,t)]R a (x'\x;t) (31) 
n 

of the social force K a (x,t) describes the reaction of subpopulation a on the behav- 
ioral distribution within subpopulation b and usually differs from K ba (x,t), which 
describes the influence of subpopulation a on subpopulation b. 
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Neglecting fluctuations, the behavior x(t) does not change if K a (x,t) vanishes. 
K a (x,t) = corresponds to an extremum of the social field V a (x,t), because it 
means 

d 

W a (x,t) = 0, i.e., —V a (x,t) = for all i G {1, . . . ,n} . (32) 

OXj, 



We will now compare our results with Lewin's (1951) social field theory. From social 
psychology it is well-known that the behavior of an individual is determined by the totality 
of environmental influences and his or her personality. Inspired by electro-magnetic field 
theory, social field theory claims that environmental influences can be considered as a 
dynamical force field, which should be mathematically representable. A temporal change 
of this field will evoke a psychical tension which, then, induces a (behavioral) compensation. 

In the following it will turn out that our model allows a fully mathematical specification 
of the field theoretical ideas, which was still to be found: 



Let us assume that an individual's objective is to behave in an optimal way with 
respect to the social field V a (x,t), that means, he or she tends to a behavior corre- 
sponding to a minimum of the social field. 

If the behavior x does not agree with a minimum of the social field V a (x,t) this 
evokes a force 

K a (x,t) = -VV a (x,t) (33) 

that is given by the gradient of the social field V a (x,t) (pulling into the direction 
of steepest descent of V a (x,t)). The force K a (x,t) plays the role of the psychical 
tension. It induces a behavioral change according to 

^*K a (x,t). (34) 

The behavioral change dx / dt drives the behavior x (t) towards a minimum x* a of the 
social field V a (x,t). When the minimum x* is reached, then 

VV a {x,t) = (35) 

holds and, therefore, K a (x, t) = 0. As a consequence, the psychical tension vanishes, 
that means, it is compensated by the previous behavioral changes. 

When the psychical tension K a (x,t) vanishes then, except for fluctuations, no be- 



havioral changes take place — in accordance with fl3~4|). The individual has reached 
an equilibrium within the social field, then. 

Note, that the social fields V a (x,t) of different subpopulations a usually have dif- 
ferent minima x*. This means that individuals of different types a of behavior will 
normally feel different psychical tensions K a (x,t). In other words, index a distin- 
guishes different personalities. 
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4 Computer simulations 



The Boltzmann-Fokker-Planck equations are able to describe a broad spectrum of 
social phenomena. In the following, some of the results shall be illustrated by computer 
simulations. We shall examine the case of A = 2 subpopulations, and a situation for which 
the interesting aspect of the individual behavior can be described by a certain position 
i on a one-dimensional continuous scale (i.e., n — 1, x = x). A concrete example for 
this situation would be the case of opinion formation. Here, conservative and progressive 
thinking individuals could be distinguished by different subpopulations. The position x 
would describe the grade of approval or disapproval with respect to a certain political 
option (for example, SDI, power stations, Golf war, etc.). 

In the one-dimensional case, the integrability conditions (^) are automatically fulfilled, 
and the social field 

X 

V a (x, t) = - [ dx' K a (x', t) - c a (t) (36) 



is well-defined. The parameter c a (t) can be chosen arbitrarily. We will take for c a (t) the 
value that shifts the absolute minimum of V a (x, t) to zero, that means, 

c a (t) := mm ( - j dx' K a (x' ,t)\ . (37) 



X(l 



• Since we will restrict the simulations to the case of imitative or avoidance processes, 
the shape of the social field V a (x, t) changes with time only due to changes of the 
probability distributions P a (x,t) (cf. (p3f)), that means, due to behavioral changes 
of the individuals (see figures [1] to [|) . 

In the one-dimensional case one can find the formal stationary solution 

Pa(x) = K(xo) ^^U^^l\ (38) 

\ XO / 

which we expect to be approached in the limit of large times t — > oo. Due to the depen- 
dence of K a (x) and Q a (x) on P a (x), equations (|38|) are only implicit equations. However, 
from (|38|) we can derive the following conclusions: 

• If the diffusion coefficients are constant (Q a (x) = Qa), (|38|) simplifies to 

P a (x) = P a (x ) exp ^~ [V a (x) + c a ] V (39) 

that means, the stationary solution P a (x) is completely determined by the social 
field V a (x). Especially, P a (x) has its maxima at the positions x*, where the social 
field V a (x) has its minima (see fig. [I]). The diffusion constant Q a regulates the width 
of the behavioral distribution P a {x). For Q a = there were no individual behavioral 
variations, and the behavioral distribution P a {x) were sharply peaked at the deepest 
minimum x* of VJx). 
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• If the diffusion coefficients Q a (x) are varying functions of the position x, the behav- 
ioral distribution P a {x) is also influenced by the concrete form of Q a (x). From (|38|) 
one expects high behavioral probabilities P a {x) where the diffusion coefficients Q a {x) 
are small (see fig. |2|, where the probability distribution P\{x) cannot be explained 
solely by the social field V\{x)). 

• Since the stationary solution P a {x) depends on both, K a (x) and Q a (x), different 
combinations of K a (x) and Q a (x) can lead to the same probability distribution 
P a {x) (see fig. |4] in the limit of large times). 

For the following simulations, we shall assume x G [1/20, 1] and use the ansatz 

e U a (x',t)-U a (x,t) 

R a (x'\x;t)= (40a) 
D a (x',x;t) 



for the readiness R a (x'\x;t) to change from x to x' (cf. ( |10a|) ). With the utility function 

"•<*■«>== "5 (Sr)'" < 40b > 

subpopulation a prefers behavior x a . L a means the indifference of subpopulation a with 
respect to variations of the position x. Moreover, we take 

Uab ^ ■= e -W-*\/r r = ^ ( 40c ) 
D a (x',x;t) ' 20 ' K J 

where R can be interpreted as measure for the range of interaction. According to (|40c|) , 
the rate of behavioral changes is the smaller the greater they are. Only small changes 
of the position (i.e., between neighboring positions) contribute with an appreciable rate. 
Figure |TJ to |5] show the respective values of R, Li, and L 2 used in the simulations. 

Note, that R and the readiness R a are different quantities. For very small values of the 
range R of interaction the diffusion coefficients Q a (x) can be neglected and the fluctuations 
play a neglible role, that means, behavioral changes are mainly given by the social field 
(see fig. 1). For greater but still small values of R, the diffusion coefficients have to 
be taken into account in order to fully understand the temporal development of the 
behavioral distribution P a (x,t) (see fig. 2). If R exceeds a certain value, the Taylor 
approximation is invalid, and the BoLTZMANN-like equations should be applied instead 
of the Boltzmann-Fokker-Planck equations. 



4.1 Sympathy and interaction frequency 

Let s a b(t) be the degree of sympathy which individuals of subpopulation a feel towards 
individuals of subpopulation b. Then, one expects the following: Whereas the frequency 
fab(t) of imitative processes will be increasing with s a b(t), the frequency fl b {t) of avoidance 
processes will be decreasing with s a b(t). This functional relationship can, for example, be 
described by 

: = £(*)«*(*). 
flit) := f a (t)(l-s ab (t)) (41) 
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with 



< s ab (t) < 1 . 



(42) 



fa(t) is a measure for the frequency of imitative processes within subpopulation a, flit) 
a measure for the frequency of avoidance processes. If we assume the sympathy between 
individuals of the same subpopulation to be be maximal, we have Suit) = 1 = ^(t)- 



4.2 Imitative processes = 1, flit) = 0) 

In the following simulations of imitative processes we assume the preferred positions to 
be xi = 6/20 and x 2 = 15/20. With 



the individuals of subpopulation a = 1 like the individuals of subpopulation a = 2, but 
not the other way round. That means, subpopulation 2 influences subpopulation 1, but 
not vice versa. 

As expected, in both behavioral distributions P a (x, t) there appears a maximum around 
the preferred behavior x a . In addition, due to imitative processes of subpopulation 1, a 
second maximum of Piix, t) develops around the preferred behavior x 2 of subpopulation 2. 
This second maximum is small, if the indifference L\ of subpopulation 1 with respect to 
variations of the position x is low (see fig. |I]). For high values of the indifference L\ even 
the majority of individuals of subpopulation 1 imitates the behavior of subpopulation 2 
(see fig. H)! One could say, the individuals of subpopulation 2 act as trendsetters. This 
phenomenon is typical for fashion. 

We shall now consider the case 



for which the subpopulations influence each other mutually with equal strengths. If the 
indifference L a with respect to changes of the position x is small in both subpopulations a, 
each probability distribution P a ix,t) has two maxima. The higher maximum is located 
around the preferred position x a . A second maximum can be found around the position 
preferred in the other subpopulation. It is the higher, the greater the indifference L a is 



However, if L a exceeds a certain value in at least one subpopulation, a socalled phase 
transition (that means, a qualitative different situation) occurs, since only one maximum 
develops in each behavioral distribution P a ix, t)\ Despite the fact that the social fields 
V a ix, t) and diffusion coefficients Q a ix,t) of the subpopulations a are different because 
of their different preferred positions x a (and different utility functions £/ a (x, t)), the be- 
havioral distributions P a (x, t) agree after some time! Especially, the maxima x* of the 
distributions P a (x, t) are located at the same position x* in both subpopulations. One 
could say, the two subpopulations made a compromise. The compromise x* is nearer to 
the position x a of the subpopulation a with the lower indifference L a (see fig. |]). 




(43) 




(44) 



(see fig. §). 
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4.3 Avoidance processes (fl(t) = 0, flit) = 1) 

For the simulation of avoidance processes we assume with x\ = 9/20 and x<± = 12/20 that 
both subpopulations nearly prefer the same behavior. Figure |5] shows the case, where the 
individuals of different subpopulations dislike each other: 

(**(*)) == ( J J ) , i-e., (£(*)) = ( J I ) - (45) 

This corresponds to a mutual influence of each subpopulation on the other. The compu- 
tational results indicate that 

• individuals avoid behaviors which are found in the other subpopulation. 

• The subpopulation a = 1 with the lower indifference L\ < L2 is distributed around 
the preferred behavior x\ and pushes away the other subpopulation! 

Despite the fact that the initial behavioral distribution P a (x, 0) agrees in both subpopula- 
tions, there is nearly no overlapping of Pi{x, t) and P2(x, t) after some time. This is typical 
of polarization phenomena in the society. Well-known examples are the development of 
ghettos or the formation of extremist groups. 

In figure |6], we assume that the individuals of subpopulation 2 like the individuals of 
subpopulation 1 and, therefore, do not react to the behaviors in subpopulation 1 with 
avoidance processes: 

(**(*)) := ( J J ) , i-e-, {f 2 ab (t)) = ( ° (46) 

As a consequence, P2ix,t) remains unchanged with time, whereas P\ix,t) drifts away 
from the preferred behavior x\ due to avoidance processes. Surprisingly, the polarization 
effect is much smaller than in figure |5|! The distributions P\(x, t) and P2ix,t) overlap 
considerably. This is, because the slope of P2ix, t) is smaller than in figure [5] (and remains 
constant). As a consequence, the probability for an individual of subpopulation 1 to meet 
a disliked individual of subpopulation 2 with the same behavior x can hardly be decreased 
by a small behavioral change. One may conclude, that polarization effects (which often 
lead to an escalation) can be reduced, if individuals do not return dislike with dislike. 



5 Empirical determination of the model parameters 

For practical purposes one has, of course, to determine the model parameters from empiri- 
cal data. Therefore, let us assume we know empirically the distribution functions P%(x, ti), 
[the interaction rates i^&(^)>] an d the effective transition rates w%(x'\x; ti) (x' ^ x) for 
a couple of times ti G {to, ■ ■ ■ ,ti}- The corresponding effective social fields V*ix,ti) and 
diffusion coefficients Q^Ax, t{) are, then, easily obtained as 

x n x 

V:(x, t{) :=-[ dx' ■ Kl(x', ti) = -T,j d < KiW, ^) ( 47a ) 

i=l 



14 



with 



K*Jx,t{) ■= jd»x , (J i -x i )wt{x , \x;t l ) 



(47b) 



and 



Q e aij (x, t{) := / d n x' (x- - Xi)(x'j - Xj)w a c {x'\x- ti) . 



(4f 



Much more difficult is the determination of the utility functions U*(x,ti), U£(x,ti), the 
distance functions D c a (x' , x; ti) , and the rates i£(t,), z/^) := ^ 6 (t/)/i 6 e (^)> *£&(*z) : = 
u ab(ti)fab(ti)- This can be done by numerical minimization of the error function 



A L 



F: =EE E i 

a=l i=0 i ; i'en 



w"(aj'|a?;t; 



1 



D a (x',x;ti) 



g a (x',x;ti) 



(49) 



for example with the method of steepest descent (cf. Forsythe et. at, 1977). In (£9), we 
have introduced the abbreviation 

g a {x\x-t{) := ^) e ^™-W0 + £ [i^(tO/?(^,*0+i4(*i)^(*.*i)]e U " ( *'' t,) - l/ ' (!BA) • 

6=1 

(50) 

It turns out (cf. Helbing, 1992c), that the rates v a (ti) have to be taken constant during 
the minimization process (e.g., u a (ti) = 1), whereas the parameters U a (x,ti), U a (x,ti), 

vlbtti) '■= Kb^i) fab(ti) and u ab(ti) '■= Kbtti) fab(ti) are to be varied. For 1/D a (x', x; t { ) one 
inserts 

1 n a (x',x;ti) 

(51a) 



D a (x',x;ti) d a (x',x;ti) 



with 



n a (x', x; t{) := w a e (x'\x\ U)g a (x' ', x; ti) [P a c (x, *,)] + w"(x\x'; ti)g a (x, x'\ ti) [P2(x f , ti) 

(51b) 

and 



d a (x',x;ti) := g a (x',x;ti)PZ{x,ti) + g a (x,x';ti)P°(x' 



(51c) 



(pl|) follows from the minimum condition for D a (x\ x;ti) (cf. Helbing, 1992c). 



Since F may have a couple of minima due to its nonlinearity, suitable start parameters 
have to be taken. Especially, the numerically determined rates u^pi) and v^b^i) have to 
be non-negative. 

If F is minimal for the parameters U a (x,ti), U a (x,ti), D a (x',x;ti), v a (ti), v^ifi) and 
Vab(ti)-> this is (as can easily be checked) also true for the scaled parameters 



U e a (x,ti) 
U^xA) 

D e a {x',x-ti) 



<{ti) 



U a {x,ti) - C a (ti) , 

U a {x,ti)-C a {ti) 

D a (x',x\tj) 

D a {ti) ' 
Vajti) 

D a (ti) ' 
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Ktit 



■*(*> - t| • (52) 



In order to obtain unique results we put 



and 



which leads to 



E^OMO^O, E^OmOesO. (53) 



^ DHx'x-t,) = S 1. ( 54 ) 

(x f ^x) (x f ^x) 



C a (t t ) := * en - , C a (U) := 2^=^ , (55) 



and 



E 



i 



Djx', x\ ti 



fl.ft) E 1 



(56) 



C a (ti) and C a (ti) are mean utilities, whereas D a (t{) is a kind of imzi o/ distance. 



The distances -D^(a?', a;; i) are suitable quantities for multidimensional scaling (Kruskal 
and Wish, 1978, Young and Hamer, 1987). They reflect the "psychical structure" (psy- 
chical topology) of individuals of subpopulation a, since they determine which behaviors 
are more or less related (compatible) (comp. to Osgood et. al, 1957). By the dependence 
on a, Dl(x', x;t) distinguishes different psychical structures resulting in different types a 
of behavior and, therefore, different "characters" (personalities). 



5.1 Evaluation of the German migration data 

It is not easy to find suitable data for the determination of the model parameters, since 
usually there only exist data for the temporal development of a certain behavioral dis- 
tribution P a (x,ti), but not for the corresponding effective transition rates w"(x'\x;ti). 
However, for a few countries all necessary data are known about migration between dif- 
ferent regions (see Weidlich and Haag, 1988). In this case, the behavior x = x means 
to live in region x G {1, 2, ... , S}, where S is the number of distinguished regions. 

In the following, the results for migration in West Germany shall be presented. West 
Germany is divided into 10 federal states and the region of West Berlin (see figure [7| and 
table [I]). The data for these S = 11 regions can be found in the Statistische Jahrbiicher 
of the years 1960 to 1985 on an annual basis (to = 1960, t\ = 1961, . . ., ti = 1985). Our 
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data analysis will assume — with respect to migration — one more or less homogeneous 
population, that means, we have A = 1, and the index a can be dropped. Figure |8| shows 
some examples for the temporal variation of the effective transition rates w e (x'\x; t{). 

Using the method of steepest descent, the minimization of the error function (^), ( |5T)| ) 
gives the following results: 

^ 2 (t;)~0, (57) 

that means, avoidance processes are negligible. The rate v e (ti) of spontaneous behavioral 
changes and the rate v le (ti) of imitative processes are depicted in figure |9|. The utility 
functions U e (x, U) for spontaneous changes and the utility functions U e (x, ti) for imitative 
processes are illustrated in figures [10] and [TT| . 

The irregulatity of the utility functions U e (x,ti) indicates that they probably fit random 
fluctuations of the migration data. Indeed, a mathematical analysis proves that the term 

U(x',ti)-U(x,ti) JJ e {x' ,ti)-U e (x,ti) 

A*)n*,u) D{x , x;ti) = s%u)n*,u) D , [x , x;ti) (58) 

only explains 5.2 percent of the variance of the effective transition rates w e (x'\x; t). There- 
fore, it makes no significant contribution to their mathematical description, and the mi- 
gration rates w e (x'\x; t) of West Germany can already be represented by the model 

w{x'\x- tt) := y{t l )—— i — = v%U)— —- — - . (59) 

D(x,x;ti) D e (x } x;ti) 



This result agrees with the model of Weidlich and Haag (1988)! Figure 12 shows 
the corresponding rate v e {ti) of spontaneous changes, and figure [13] depicts the utility 
functions U e (x,ti). The distances D e (x' ,x;ti) can be calculated from the formulas (|5"T|), 
(|52|), and (p6[). They are not only a measure for the mean geographical distances, but 
also for transaction costs (e.g., removal costs) and psychical differences (of the language, 
mentality, etc.). 

The replacement of the time dependent distances D e (x',x;ti) with the time independent 
values Dl(x',x) defined by 

DZ(x',x) '~ T+if^D^x^x-ti) (60) 

(see table ^) allows for a further model reduction. The optimal value of the rate of spon- 
taneous behavioral changes is, then, given by 



x,x ■■ 



Although the reduced model 



e U°(x',ti)-U°{x,ti) 



Hx'M := » e M D%{x ^ x) (62) 

only needs (S + 1) • (L + 1) + S ■ (S — l)/2 = 367 variables for the description of S ■ (S — 1) ■ 
[L + 1) = 2860 effective transition rates w e (x'\x; ti), it attains a very good correlation of 
0.984 with the empirical data. For a more detailled discussion of this model, see Weidlich 
and Haag (1988). 
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6 Summary and outlook 



In this article, a behavioral model has been proposed that incorporates in a consistent way 
many models of social theory: the diffusion models, the multinomial logit model, Lewin's 
field theory, the logistic equation, the gravity model, the Weidlich-Haag model, and the 
game dynamical equations. This very general model opens new perspectives concerning 
a theoretical description and understanding of behavioral changes, since it is formulated 
fully mathematically. It takes into account spontaneous (or externally induced) behavioral 
changes and behavioral changes due to pair interactions. Two important kinds of pair 
interactions have been distinguished: imitative processes and avoidance processes. The 
model turns out to be suitable for computational simulations, but it can also be applied 
to concrete empirical data. 



6.1 Memory effects 



The formulation of the model in the previous sections has neglected memory effects that 
may also influence behavioral changes. However, memory effects can be easily included 
by generalizing the BOLTZMANN-like equations to 



d r 

-Pa(x,t) = dt'Y, [wU'(x\x'-,t')Pa(x',t') - W?_ t ,( X '\x-,t')P a (x,t') 



to ^ 

with the effective transition rates 



o := (x'\x-, O + EE E «&* y>> y> f - 

b=iyeny'en 

and generalizing the Boltzmann-Fokker-Planck equations to 

t l ~ l 

1 " d d 



+^E 



gg(aj,OPa(*,0 



2 i%_ dx t dxj 

with the effective drift coefficients 

K£*(xJ) := jd^'ix'.-x^wUiAx-t'), 
n 

the effective diffusion coefficients 

Q*7/(*>0 : = J^x'M-x^-x^wU'ix'WJ), 



and 



wU{x'\x-t') :=w t - t '(x > \x;t') + 't I d n y I d n y' w'J \x' ,y>\x,y;t>)P b (y,t') . (64d) 

6=1 J 



(63a) 



(63b) 



(64a) 



(64b) 



(64c) 



n n 
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Obviously, in these formulas there only appears an additional integration over past times 
if (Helbing, 1992c). The influence of the past results in a dependence of wf_ t ,(x'\x;t'), 
K^*' (x,t f ), and Q t ^f(x,1f) on (t - t'). The BoLTZMANN-like equations @ resp. the 
Boltzmann-Fokker-Planck equations used in the previous sections result from 
(|63|) resp. fl64[) in the limit of short memory. 



6.2 Group dynamics 

The force model described in section |] can serve as a new mathematical modelling concept. 
For example, it has successfully been applied to the simulation of pedestrian behavior (cf. 
Helbing, 1991). 

Another interesting field is an application of the force model to group dynamics and group 
formation. In this case we take A = N, that means, each subpopulation consists of one 
individual only (N a = 1). Since N a 3> 1 is violated, then, the temporal change dP a (x, t)/dt 
of P a (x, t), which describes the probability of individual a to show the behavior x at time t, 
is additionally subject to fluctuations (cf. Helbing, 1992a, c). The interaction rates = 
v a b are related to the adjacency matrix and describe the social (interpersonal) network 
(Burt, 1982). Moreover, the sympathy matrix (s ab (t)^j is affected by the social processes. 
The crucial task for simulating group dynamics is, therefore, to set up equations for the 
temporal change of s a b(t). The topic of group dynamics will be treated in a forthcoming 
paper. 
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Figure 1: Imitative processes in the case of one-sided sympathy and low indifference 
L a with respect to behavioral changes. 




Figure 2: As figure 1, but for high indifference L a with respect to behavioral changes. 




Figure 3: Imitative processes for mutual sympathy and low indifference L a in both 
subpopulations. 




Figure 4: As figure 3, but for high indifference L 2 in subpopulation 2. 




Figure 5: Avoidance processes for mutual dislike of both subpopulations. 



Figure 6: Avoidance processes for one-sided dislike. 




Figure 7: The subdivision of West Germany into eleven federal states (from Weidlich 
and Haag, 1988). 
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Schleswig-Holstein 
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Niedersachsen 
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Bremen 


O 
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Nordrhein-Westfalen 
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Hessen 
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Rheinland-Pfalz 
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Baden-Wurttemberg 
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Bayern 
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10 


Saarland 


/ 


11 


Berlin 



Table 1: The eleven federal states of West 
in the following figures. 
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Figure 8: Time dependence of some effective transition rates of migration in West Ger- 
many. 
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Figure 9: Rate v e (ti) of spontaneous behavioral changes (□) and rate z/ le (t;) of imitative 
processes (A) calculated from the migration data of West Germany. 
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Figure 10: Utility functions for spontaneous removals from one federal state to another in 
West Germany. 
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Figure 11: Utility functions for removals from one federal state to another due to imitative 
processes. 
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Figure 12: Rate of spontaneous migration in West Germany for the model of Weidlich 
and Ha AG. 
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Figure 13: Utility functions of the federal states of West Germany for the model of Wei- 
dlich and Haag. Note, that the utility of West Berlin (/) is extremely dependent on 
changes of the political situation. For example, there is an remarkable increase of the 
utility after the erection of the Berlin wall in the year 1961. 



X 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


1 


- 


0.16 


0.56 


1.19 


0.99 


1.93 


3.01 


1.79 


2.13 


8.35 


1.26 


2 


0.16 


- 


0.44 


1.65 


1.58 


2.18 


5.03 


2.34 


2.61 


12.89 


1.55 


3 


0.56 


0.44 


- 


0.25 


0.44 


0.90 


2.13 


1.26 


1.55 


6.38 


0.79 


4 


1.19 


1.65 


0.25 


- 


1.91 


3.10 


6.16 


3.31 


4.30 


17.37 


2.58 


5 


0.99 


1.58 


0.44 


1.91 


- 


0.65 


0.53 


0.79 


0.93 


2.56 


0.92 


6 


1.93 


2.18 


0.90 


3.10 


0.65 




0.47 


0.61 


0.75 


2.23 


1.15 


7 


3.01 


5.03 


2.13 


6.16 


0.53 


0.47 




0.59 


1.39 


0.54 


2.33 


8 


1.79 


2.34 


1.26 


3.31 


0.79 


0.61 


0.59 




0.38 


1.56 


1.14 


9 


2.13 


2.61 


1.55 


4.30 


0.93 


0.75 


1.39 


0.38 




3.49 


1.09 


10 


8.35 


12.89 


6.38 


17.37 


2.56 


2.23 


0.54 


1.56 


3.49 




5.07 


11 


1.26 


1.55 


0.79 


2.58 


0.92 


1.15 


2.33 


1.14 


1.09 


5.07 





Table 2: Time independent effective distances Dl(x',x) between the eleven federal states 
of West Germany. 



